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Abstract. We perform calculations of pre-main sequence evolution of stars from 1 to 85 Mq with growing accretion 
rates M. The v alues of M are taken equa l to a constant fraction / of the rates of the mass outflows observed 



by Churchwell (1998) and Henning (2000). The evolution of the various stellar parameters is given, as well as 



the evolution of the disc luminosity; electronic tables are provided as a supplement to the articles. Typically, the 
duration of the accretion phase of massive stars is ~ 3 • 10 5 yr. and there is less than 10% difference in the time 
necessary to form a 8 or 80 Mq star. If in a young cluster all the proto-stellar cores start to accrete at the same 
time, we then have a relation M(t) between the masses of the new stars and the time t of their appearance. 
Since we also know the distribution of stellar masses at the end of star formation (IMF), we can derive the star 
formation history N(t). Interestingly enough, the current IMF implies two peaks of star formation: low mass stars 
form first and high mass star form later. 
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1. Introduction 

Two scenarios have been proposed to explain the forma- 
tion of massiv e star s. Star coalescence (Bonnell et al. 1998 ; 
Stahler et al . 1998 ) and the accret ion sc enario (Beech & 
Mitalas 1994 ; Bernasconi & Maeder 1996 ). The merging of 
stars requires very high stella r dens ities (> 1 4 sta rs pc~ 3 ) 
to be efficient (Bonnell et al. |1998[ Henning p00l| ) . 

The imaging of new-born stars with the HST (Burrows 
et al. [1996D , the 1RAM interferometer (Dut rey |1998| ), the 
VLA (Wilner & Ho fc Kastner & Rodriguez ^OOOp and the 
VLT (Brandner pOOOj ), among others) clearly shows discs 
(or remnants of discs) , favouring the accretion scenario for 
low and intermediate mass stars. The jets and/or outflows 
that are sometime observed are closely related to accretion 
discs; they can contribute to solve the angular momentum 



problem (Tomisaka 200C). 

Accretion rates of the order of 10 -5 to 10~ 4 MQyr -1 
explain very well low mass stars (Palla 1998). But such 



values fail to describe the formation of massive stars, be- 
cause with such rates the massive stars would leave the 
ZAMS (i.e. would have burned a significant part of their 
hydrogen ) befo re being fully formed, as pointe d out by 
Nakan o ( |l988| ) and Norberg & Maeder (|2000| ). Nakano 
( 200C ) studied the implications of a huge constant accre- 
tion rate of » 10 -2 Mq yr -1 for an observed 10 5 L( 
tostar. 



The reality could well lie between these two models. 
Can an accretion rate that depends on the stellar mass 
(or luminosity) describe the star formation from low to 
very high masses ? This is the question of the growing 
accretion rates. 

A nice observational relation between the outflow mass 
rates and the stellar bolometric luminosities in ultra com- 



pact Hn regions was established by Churchwell (1998) 
over the ranges 1 to 10 6 Lq and 10~ 6 to 10~ 2 Mq yr 



Independent observations by Henning et al. (2000) confirm 
its validity. 

Churchwell's relation can bring some new constraints 
for the scenarios of the formation of massive stars. A cru- 
cial question is to know if the accretion scenario is com- 
patible or not with this relation. The relation can bring 
something new to the description of the pre-main sequence 
evolution of stars and for the explanation of the rapid birth 
of massive stars. The relation can also be a used as a strong 
check for the computer simulations of discs around stars 
and the jets/outflows, and offers the possibility to fix un- 
certain parameters of these models (e.g. viscosity). 

The principal properties of stars with accretion discs 
on the birthline are discussed. Tracks for stars leaving the 
birthline are calculated. Links between the IMF and the 
processes that end the stellar mass evolution are also pre- 
sented. 
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Fig. 1. The logT c ff - logL diagram showing birthlines for 3 values of / and for the model of N&M. The ZAMS is 
taken from Schaller et al. (1992). Observations of protostars from the compilation by N&M are shown with discs 
(respectively circles) if M is more (respectively less) than 10~ 5 MQyr _1 ; triangles mean unspecified rates. 



2. The simulations 

The matter which escapes from the disc (in the central 
region) is divided in two major flows: one part falls on the 
star while the other part is ejected in the form of jets and 
outflows. Let us denote by M^isc the mass rate transferred 
inwards from the disc, by M the mass rate accreted by the 
star, and by M out the mass rate in the jets/outflow. 

The relation, found b y Ch urchwcll (|l99S|) and con- 
firmed by Henning et al. ( 2000| ), between the bolometric 
luminosity L of stars and the outflow mass rate M out is 
single valued (at least in the error bars), monotonously 
increasing and valid over a very wide range of luminosity: 



from 1 to 10 6 L(7). A polynomial fit of the observations by 



Churchwell (199S) gives 



log M ° ut - = -5.28 + log A • (0.752 - 0.0278 log A). 



Unfortunately, 
available between 
ratio between the 
some suppositions 
able to model the 
M = /M out , M = 
that /=//(!-/)• 



there is at present no similar relation 
L and M; one must parameterize the 
accretion and outflow fluxes and do 
about their proportion before being 
evolution of the central star. Let be 
/ Afdisc and M out = (1 - f)M disc so 
A variety of physical parameters could 
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have effects on / (or /): chemical compositions, metallici- 
ties, magnetic fields, rotational velocity of the star. Some 
other parameters have certainly a role to play, but are not 
yet firmly defined, like the dust opacity. Therefore, obtain- 
ing a theoretical value of / is very difficult. Nonetheless, 
Tomisaka ( |1998| ) found / « 1/3 with his numerical MHD- 
simulation of a low mass class protostar. Shu et al. (1998) 
deduced a similar value with their model of X-wind mag- 
netic configuration, also for a low mass star. From his 
measurements of the far-infrared luminosity of a B star 



and the mass of its outflows, Churchwell (1998) estimated 
a slightly smaller value, of the order of / « 0.15. 

Due to the great number of changing physical param- 
eters which can have an effect on /, / is probably non- 
constant during the evolution of accreting stars. The last 
three values of / cited above suggest that / could be a de- 
creasing function of the stellar mass, but this is far from 
being firmly established. For this preliminary study, we 
simplify this huge astrophysical problem, retaining only a 
ratio / that is constant during the whole accretion period. 

The " best" / is presently unknown. A way to assign a 
more or less realistic value to / is as follows. The birthline 
is the path in the HR diagram that continuously accret- 
ing stars follow. Stars on the birthline arc difficult to ob- 
serve, because they are hidden by the accretion disc and 
a dense cocoon of interstellar matter. The star becomes 
visible once a sufficient amount of the surrounding mat- 
ter has been dispersed. It follows that the young stars are 
observed generally between the birthline and the zero age 
main-sequence (ZAMS). Before hydrogen ignition, gravi- 
tational contraction and deuterium burning are the main 
energy sources for the star; the localization of the birthline 
depends on the deuterium accretion rate. We will choose 
/ as the value for which the computed birthline fits "at 
the best" the upper envelope of the observations of young 
stars in the pre-main s equen ce part of the HR diagram (as 
in Norberg & Maeder |200C| , denoted hereafter N&M). 

The evolution of the pre-main sequence models with 
growing accretion r ate is computed with the Geneva code 
(Meynet & Maeder 2000] ), under the assumption of a non- 
rotating star, at Z — 0.02 and at deuterium concentration 
5 • 10~ D . Winds are not taken into account. 

The computations begin with a fully convective 
0.7 Mq protostar which is taken to be 7 • 10 5 yr old, 
in agreement with models of low mass star formation 
(Palla 1998| ). The mass of the protostar grows at a rate 
M = max(10- 5 M yr" 1 , f M mt ). 

We compare the models with accretion rates param- 
eterized by the bolometric luminosity of the star to the 
model of mass parameterized accretion rate found by 
N&M: M = 10~ 5 -max(l, A/) 15 where units are Mq and 
Mq yr _1 . The M- and L-parameterized accretion rates 
are different in the sense that the first one depends only 
on one parameter (the mass) of the star, and the sec- 
ond one depends on many involved quantities (mass, pro- 
files of chemical composition and internal energy, etc.) and 
thus is sensitive to the history of the star (e.g. opacity via 
chemical abundances, diffusive processes). For this rea- 
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Fig. 2. A birthline in logT c ff - logi. Some key phases 
are highlighted: 0: initial conditions - 1: central ignition 
of deuterium - 2: central gravitation dominates over 2 H- 
burning - 3: start of shell deuterium-burning - 4: end of 
large scale convection - 5: end of thin central convection - 
6: central 2 H -exhaustion - 7: central ignition of hydrogen 
and restart of core convection - 8: arrival on the ZAMS - 
9: star starts evolving away from the ZAMS. 



son, the memory of the initial luminosity of the star in 
the M-parameterization is rapidly lost. This is not the 
case for the L-parameterized accretion; if the initial star 
is over-/under-luminous, the luminosity of the accreted 
deuterium is also above/below the normal value and the 
birthline is also over- /under- luminous. This sensitivity to 
the initial conditions for the L-parameterized accretion re- 
mains small with respect to the possible variations of the 
birthlincs with /, for stars initially on the Hayashi band. 
Processes that depend on their past can " easily" os cillate 
or even be chaotic, for physical reasons (Hairer et al. 1987) 
or for numerical reasons (Berge et al. |1988 ). Apparently 
the calculations do not reveal strong signatures of insta- 
bilities, except maybe some oscillations near the surface 
of accreting stars above 22 Mq. 



3. Results and discussions 

The first step consists of the estimate of a plausible /. The 
second step is the prediction of physical parameters for the 
stars and discs using the / found. Connections with the 
IMF is postponed to the next section. 

Before hydrogen ignition, the main energy sources for 
the star are gravitational contraction and deuterium burn- 
ing. The localization of the birthline is then governed es- 
sentially by the deuterium accretion rate on the star. The 
higher M, the higher the birthline in the HR diagram. 
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Table 1. Properties of stars on the birthline with / = 
0.5: age, luminosity, mass, effective temperature, radius, 
relative upper radius of the convective zone. The ratio of 



disc and 


star luminosities is also 


given. 






t/yr 
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Fig. 3. Location of the convective and radiative zones. 
The relative radius of the boundaries between the convec- 
tive and radiative zones are drawn for the birthline with 
/ = 0.5. For the other birthlines, only the beginning of 
the fully radiative state is shown. The small labels denote 
the evolutionary stages shown in Fig. ||. 



Figure [l] shows birthlines for different values of /, and 
the observations compiled in N&M. While / = 0.3 seems 
to be a little too low, / £ {0.4, 0.5} give reasonable en- 
velopes to the observations. Two statistical biases in the 
estimate of the envelope are possible, with consequences 
for the determination of /. 1) The stars which constrain 
the birthline arc usually all on one side of the birthline (a 
nonsignificant exception is nevertheless visible in Fig. |J). 
Because the density of stars near the birthline is small, the 
probability that one misses a crucial observation is rela- 
tively large. The effect of this bias is an underestimated 
luminosity of the envelope of the observation and thus 
also of the birthline. Unfortunately, this occurs in the up- 
per pre-main sequence part of the HR diagram where the 
sensibility of the computed birthline to / is the greatest. 
The deficiency of stars in this region is due to the very 
short time spent by a star in this region, as will be shown 
in Fig. ||. 2) The stars leave the birthline faster than they 
reach the ZAMS. In other words, the median of the po- 
sitions of the stars for which the error bars constrain the 
birthline (by their envelope) is shifted towards the ZAMS; 
this effect vanishes as error bars go to zero. With these 
two facts in mind, our preference is for the highest value: 
/ = 0.5, i.e. / = 1/3. 

The ratio / = A/out /Miisc supported by Fig. |l| is in 
good agreement with theoretical estimates of / = 1/3, as 
mentioned above. 



Properties of stars on the birthline and the relative 
luminosity of the disc to that of the star, for / = 0.5, 
are summarized in Table |], while Fig. |^ locates different 
phases of the star evolution in the HR diagram. 

The 2 H -burning begins at the center of the star. 
Nascent stars of mass below 3.8 Mq (log(M/M Q ) = 0.58) 
are fully convective (Fig. |J), and the deuterium accreted 
is carried inside the star where it is burnt. As the temper- 
ature rises, the deuterium burns farther from the center. 
At 3.8 Mq, a radiative core is created, and the accreted 
deuterium is then burnt in a shell at the periphery of this 
core; the star inflates (second hump in Fig. ||) and the 
luminosity grows. If the relation of the L-parameterized 
accretion is right, and in particular its single valued char- 
acter, the shell 2 H -burning could play a major role in 
the beginning of star formation, raising the stellar lu- 
minosity and inducing via Churchwell's relation, a real 
boost in the accretion rate, shortening the time needed to 
reach high mass stars. The convective 2 H -burning shell 
becomes thin at the surface of the star around 5.6 Mq 
(log(M/MQ) = 0.75) (slightly depending on the mass ac- 
cretion rate); the star becomes fully radiative. 

The luminosity of the 2 H -burning is governed by the 
accretion rate; the star inflates to balance a excess of 
power, or contracts during a shortage of deuterium. The 
R(M) relation (Fig. |) is thus very sensitive to M during 
the pre-main sequence evolution of the star and especially 
during the fully radiative phase. 
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Fig. 5. Mass evolution for stars on the birthlines. 

The consequences of the luminosity's growth during 
the shell 2 H -burning are clearly visible in the logt - 
logM diagram (Fig. |5|) in the form of an elbow near 
log(M/Af ) = 0.8 (6M ) with a rather abrupt change 
of slope in the curves, for the L-parameterized accre- 
tion. Below the bend, the much smoother curve of the 
M-parameterized accretion of N&M is a little bit higher 
than the one for the ^-parameterized accretion, and lower 
afterwards. Rather than physical, the origin of this dif- 
ference resides in the choice of the parameterization for 
M. Because both laws were adjusted on the same obser- 
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Fig. 7. Relative luminosity of disc to that of the star as 
a function of stellar bolometric luminosity for stars evolv- 
ing on the birthlines. For the model of N&M, / = 0.5 is 
assumed. 

vations, they globally look alike and make similar predic- 
tions, but they can locally differ. 

One half of the gravitational energy gained during the 
rotational infall is transformed into kinetic energy; the 
other half is used to heat up the disc and is eventually ra- 
diated. Chemical energies (sublimation of dust and grains, 
ionization) are small with regards to gravitational energy 
and are neglected in the computations. The disc luminos- 
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Fig. 8. Tracks in logT c ff - logL followed by stars that leave the birthline with / = 0.5. Couples of numbers represent 
the age i s t op of the star (in units of 10 3 yr) at the stop of the accretion, and its mass (in Mq) on the ZAMS. For the 
stars which are no longer on the birthline, the small dots at logT c ff ~ 3.75 correspond to the beginning of the fully 
radiative interiors; the small dots close to the ZAMS correspond to the appearance of the convective cores. 



ity, Ldisc, is thus of the order of G M M disc / (2 R), as given 



by Hartmann (1998a, §5.4). It appears in Fig. Qthat dur- 
ing the accretion phase, the luminosity of the disc is always 
greater than the star luminosity. On the birthline, the star 
encounters three phases of nearly constant luminosity. In 
these phases, the radius of the star shrinks while the mass 
continues to grow; according to the previous formula, the 
luminosity of the disc grows. That explains the three ver- 
tical steps in the logL - log(L/Ldisc) diagram. The step 
around log(L/-Lo) = 4.8 corresponds to the arrival of the 
star on the ZAMS and is the only feature really sensitive 
to /. The probability of having good observations of stars 



in this region of the diagram is very small and we conclude 
that the relative luminosity of the disc to that of the star 
cannot help us to fix /. The smaller variability of the ratio 
of disc to star luminosities for the L-parameterized accre- 
tion than for the M -parameterized accretion has the same 
origin as in the log t - log M diagram. 

After the accreting star has joined the ZAMS, at 
22 Mq, it evolves in a very similar manner to the one cal- 
culated in N&M. They used M oc M 1 ' 5 and the linear fit in 
Churchwell's data implies M oc M 4 using the L(M) re- 
lation on the ZAMS. An accretion rate of 8- 1CT 3 Mq yr" 1 
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is obtained for a 85 Mq star. Though it seems very high, 
this rate is still in the permitted domain of accretion with 
spherical symmetry (Wolfire & Cassinelli 1987); if the ac- 
cretion is made via a disc, even higher rates are permitted 
due to the less important cross-section of infalling mate- 
rial to radiations from the star and from the shock, and 
to the enhanced cooling (Nakano 198S; Henning [2001 ). 

Once the disc is no longer fed by its environment, for 
example by the disappearance of the parent cloud by the 
outflows of the star or jets/ouflows and the r adiati on of 
neighbouring hot luminous stars (Nakano et al. 1995 ), the 
disc get thinner by accretion and evaporation. The tran- 
sition from active to passive discs is probably a smooth 
function of time. As its mass decreases, the residual disc 
plays a smaller and smaller role in the evolution of the 
star. Both the mass of discs and the typical duration of 
the transition are still very uncertain; in the computations, 
we neglect their effects, using an instantaneous stopping 
of the accretion at time t stop . Figure || shows the result- 
ing theoretical tracks for some values of t s top (m units of 
10 3 yr). 

The following considerations can partly justify the as- 
sumption that from a certain stage, the disc has no more 
significant influence on the central star. For low and in- 
termediate mass stars, the disc is usually lighter than 
the central star (Hartmann [19984 §§ 6 - 3 and 9 - 7 > |1998b ; 
Bcrnasconi 1997). Observations of the residual mass trans- 
fer rate of the disc around T Tauri stars is described 
by Hartmann ( 1998b| ); they can be approximately repre- 
sented by logM£g S c pa 0.4 - 1.4 logf (units are Mq yr -1 
and yr), where t is a relative scale for the ages. Integrating 
M|i| c from the i corresponding to a disc mass transfer rate 
(1 + /) M out at i st0 p to infinity, one finds a crude estimate 
of the mass of the disc at the moment of the stopping 
of the accretion: Mdi sc ~ 5 Mdisc ' 3 (units are Mq and 
Mq yr -1 ). Mdisc is probably out of the validity domain 
for M" s c , but /M di! 



remains less than 5% of the mass 



of the star, for M > 2 Mq. Bernasconi ( |1997[ ) studied 



the residual mass accretion of dying discs, but his models 
contain many arbitrary parameters; for intermediate mass 
stars, the difference between the tracks with and without 
residual accretion are nonetheless much smaller than the 
error bars of the present day observations. 

The time t st0 p corresponds approximately to the time 
at which the cocoon of interstellar matter vanishes, thus 
the star becomes visible in the optical bands only after 
this moment; the location of the first optical visibility of 
stars above 22 Mq is very close to the ZAMS. 

For the stars studied in this paper (mass on the ZAMS 
> 1.19 Mq), there are four regions in the pre-main se- 
quence part of the HR diagram, as shown in Fig. |8[ 1) 
For logr off below « 3.68 (i.e. 4.8 • 10 3 K), the stars are 
fully convective; they leave the birthline with a vertical 
Hayashi-like displacement in the HR diagram. No more 
accreting stars between logT ff sa 3.68 and logT c ff sa 3.76 
(i.e. 5.8 • 10 3 K) have a radiative core surrounded by a con- 
vective 2 H -burning shell; if the star already has a radiative 



core before the end of accretion, the track closely follows 
the birthline (with an increasing luminosity) for a short 
period of time after t s top; if a radiative core is absent, one 
is swiftly created after the end of accretion and the star 
leaves the birthline with a decreasing luminosity. 3) For 
log T c ff above f=a 3.76, the stars are fully radiative, except 
on 4) a small band above the ZAMS where stars already 
have a convective H -burning core. 

Electronic tables for the tracks are available via the 
anonymous ftp server of the Geneva observatory; the file 
is located in directory 



ftp : //obsftp.unige.ch/pub/evol/prems/AA2001BM 



4. Relation with the IMF 

The initial mass function (IMF) represents the density 
of probability that a star ends its accretion with a given 
mass. The IMF results from many physical phenomena 
(fragmentation of the parent cloud, vanishing of the reser- 
voir of interstellar matter, accretion rate and its efficiency, 
binarity, etc.). But under the assumption that the accre- 
tion rate does not depend on the local conditions of the 
cloud and disc, the IMF can be seen as the result of the 
mechanisms that end the feeding of the disc with inter- 
stellar matter (in the case of a small ratio of disc to star 
masses) . 

The IMF is usually noted £(m). £,(m) dm is the relative 
number of new-born stars in the mass interval dm centered 
on m; the IMF is normalized so that J M max £(m)dm = 
1. The assumptions made for the models of this paper 
imply that the star mass evolution on the birthline is a 
single valued and monotonously increasing function M(t). 
Thus, from the observed IMF and using the M(t) relation 
(as in Fig. |]) , one can compute the density of probability 
-P (*stop) of the age i s top of stars when the accretion is 
abruptly stopped: 

P(t s to P ) = (£(M(i))dM(i)/dt)| t=tBtop . 

The result is shown in Fig. [| using the IMF of Kroupa 
( [2001P : £(m) = /3m" Q where 



(0.30; 1.99) , if m G [M min = 0.01, 0.08), 
(a; (3) = { (1.30; 0.159) , if m £ [0.08, 0.50), 

(2.30; 0.0795) , if m £ [0.50, M max = 100.] 



The IMF of Scalo (JT99|) differs from Kroupa's one prin- 
cipally in the < 0.1 Mq region; both IMFs give similar 
results. 

For the L-parameterized accretion, the break in the 
mass evolution graph located near 6 Mq is great enough 
to generate a sharp hump (an order of magnitude) 
in the density of probability of t s top, f (istop)- For 
the M-parameterized accretion of N&M which gives a 
monotonously decreasing P(i s top), a bump is also distin- 
guishable, but it is smoothed by the nature of the law used. 
In comparison to the L-parameterized accretion rate, the 
M-parameterized accretion rate is indeed a little higher 
for masses in the range 1.6 to 6Mq and lower for the 
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Fig. 9. Density of probability of t stop for mass accre- 
tion rates with three values of /, and the N&M one. For 
istop < 7 ■ 10 4 yr [M < 0.7 Mq), an accretion rate of 
10~ 5 Mq yr _1 is assumed. 

other stars (Fig. |J). The presence of the peak for the 
L-paramctcrizcd accretion docs not depend on the ratio 
/ = M / M out of accreted and ejected masses, if this quan- 
tity does not evolve with the mass and luminosity of the 
star. If theses peaks are confirmed, it would be most inter- 
esting to investigate the various physical processes at their 
origin. One possible way to explain theses curves could be 
as follows. For stars below 6Mq, the fractal parent cloud 
is probably dissipated mostly by neighbouring stars. For 
stars above Q Mq, the star and disc luminosities could 
grow high enough to scatter the parent cloud. 

5. Conclusions 

Numerical experiments have been done with mass accre- 
tion rates proportional to the luminosity parameterized 
outflow rate 



log- 



= -5.28 



log-— 

Lq 



(0.752 - 0.0278 log— ) 

Lq 



by Churchwell (1998) and confirmed by 



M yr- 
established 

Hcnning et al. (200C). Until complete theoretical values 
of the ratio of accreted on the star and ejected masses, 
/ = M/M out , are computed taking into account stellar 
rotation and magnetic field, chemistry, radiative transfers, 
etc., / can be seen as a unknown parameter. 

1. The comparison of theoretical birthlines from differ- 
ent / to observations of young star constraints /. For 
these experiments, we chose constant values for / and 
our favorite choice is / = 0.5. This value permits a 
reasonable fit of the observations of accreting stars. 



Only a fraction / = //(l + /) = 1/3 of the mass 
transferred from the disc in the central region accretes 
on the star. This ratio is in relative agreement with 
numerical predictions by Tomisaka ( 1998| ) (1/3) and 
the theoretical model of Shu et al. (1998) (also 1/3), 
both cases are for low mass stars. Observations by 



Churchwell (1998) suggest a smaller value (of the or- 
der of 0.15) for a B star; this is not seen as a strong 
disagreement. 

The relative lack of observations in the upper pre-main 
sequence part of the HR diagram is due to the very 
fast evolution of the stars just before 1 H -ignition and 
their arrival on the ZAMS; there are few chances to 
better constrain / using the HR diagram. Apparently, 
the ratio of disc to star luminosities cannot help to fix 
/• 

The R{M) relation is very sensitive to the mass ac- 
cretion rate during the 2 H -burning; the biggest dif- 
ferences are located in the range of 7 to 15Mq, i.e. 
when the star is fully radiative and the 2 H burns near 
the surface; for a given mass, the greater the rate, the 
bigger the radius. 

The M-parameterized accretion rate of Norberg & 
Maeder (200C) gives essentially the same birthline in 
the HR diagram as the L-parameterized accretion. The 
mass evolution of stars with accretion rate parameter- 
ized by M is smooth in logt - logM; on the contrary, 
the evolution for the mass with the L-parameterized 
accretion shows a sharp transition near 6 Mq . This 
is related to the increasing radius and luminosity of 
the star during the shell 2 H -burning. This dramatic 
change, as seen in the diagram of Fig. ^ has cer tainly 
a counterpart in the shape of the IMF. Kroupa ( 2001 ) 
guesses that the slope of the observed IMF could have 
a break near 6 Mq, just in the elbow of M(t). 
At 85 Mq, the accretion rate is about 8 TO -3 Mq yr _1 . 
This tremendous rate is not forbidden by star and 
shock luminosities, even in the unfavorable spherical 
case studied by Wolfire & Cassinelli fll987| ). 
The time needed to form a 85 Mq star is about one 
quarter of Myr for / = 0.5 and one third of Myr for / = 
0.3. In all cases, it is much shorter than the hydrogen 
burning timescale. 



Growing accretion seems to be a valuable scenario to 
explain the formation of massive stars, but the field for 
theories about the physical backgrounds of the accretion 
processes is still widely open. There are open questions for 
which this prospective work could not give answers. Is the 
hypothesis of the constancy of / valid for all stars on the 
birthline ? How can one explain such a high M ? Is M(L) 
really single valued, as suggested by the observations ? 
What occurs in the accreting stars observed far from the 
birthline (for 5 < log(L/Z/Q) < 6) ? Observational, the- 
oretical and numerical works should be done in order to 
better understand the mechanisms of the huge accretion 
processes. 
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